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First theoretical and numerical results on the global structure of the energy shell, the Green function 
spectra and the eigenfunctions, both localized and ergodic, in a generic conservative quantum 
system are presented. In case of quantum localization the eigenfunctions are shown to be typically 
narrow and solid, with centers randomly scattered within the semicircle energy shell while the 
Green function spectral density (local spectral density of states) is extended over the whole shell, 
but sparse. 
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One of the main results in the study of the so-called 
quantum chaos has been the discovery of quantum dy- 
namical localization as a mesoscopic quasi-classical phe- 
nomenon 0]. This phenomenon has been widely studied 
and confirmed by many researchers for dynamical models 
described by maps. Contrary to a common belief, maps 
describe not only time-dependent systems, but also con- 
servative ones (in the form of Poincare' maps). On the 
other hand, to our knowledge, there are no direct studies 
of quantum dynamical localization in bounded conser- 
vative models; moreover, the appearance of dynamical 
localization in such systems due to quantum effects is 
challenged by some researchers. The existence of local- 
ization in conservative systems would restrict quantum 
distributions to smaller regions of phase space than clas- 
sically allowed, and would therefore introduce significant 
deviations from ergodicity. 

We have addressed this problem on the Wigner Band 
Random Matrix (WBRM) model, which was introduced 
by Wigner 40 years ago |^ for the description of complex, 
conservative quantum systems like atomic nuclei. Due to 
severe mathematical difficulties, the random matrix the- 
ory (RMT) immediately turned to the much simpler case 
of statistically homogeneous (full) matrices, for which im- 
pressive theoretical results have been achieved (see, e.g., 
Refs. [^). However, full matrices describe local chaotic 
structures only, and this limitation is often inacceptable, 
for instance in the case of atoms . 

Generally speaking, RMT is a statistical theory of 
systems with discrete energy (and frequency) spectrum. 
Since the latter is a typical property of quantum dynam- 
ical chaos Q, RMT provides a statistical description of 
quantum chaos and, what is very important, one which 
does not involve any coupling to a thermal bath, which 
is a standard element in most statistical theories. More- 
over, a single matrix from a given statistical ensemble 
represents the typical (generic) dynamical system of a 
given class, characterized by a few matrix parameters. 
This makes an important bridge between dynamical and 



statistical description of quantum chaos. 

To the extent that Band Random Matrices can be 
taken as models for generic few-freedoms conservative 
systems which are classically strongly chaotic (in partic- 
ular ergodic) on a compact energy surface, the results 
presented in this Letter provide the first characterization 
of the properties of quantum chaos in momentum space 
for quantum systems of this class. 

We consider real Hamiltonian matrices of a rather gen- 
eral type (more specific random matrix models have been 
recently proposed in j7j) 



Hr, 



(m, n ~ 1, .., N) 



(1) 



where off-diagonal matrix elements Vmn — Vnm are sta- 
tistically independent, Gaussian random variables, with 
< v„in >~ and < u^, if |m — n| < 6, and are 

zero otherwise. In a classical picture, WBRMs like (|l|) 
would correspond to classical Hamiltonians of the form: 



H = Ho + V 



(2) 



where the perturbation V is usually assumed to be suffi- 
ciently small, while the unperturbed Hamiltonian Hq is 
completely integrable. In the quantum model the matrix 
(^ is given in the basis of the unperturbed eigenstates 
of Hq. Correspondingly, the fluctuations of unperturbed 
energy levels e„ are taken as Poissonian. Although in 
completely integrable quantum system there is a quan- 
tum number for each freedom, we suppose that the un- 
perturbed states are ordered according to increasing en- 
ergy, and we thereby label them by a single number n. 
The most important characteristic of WBRM is the av- 
erage level density p: 



(3) 



Here and below, the averaging is understood either over 
disorder (that is, over many random matrices) or within 
a single, sufficiently large, matrix. Both ways are equiva- 
lent owing to assumed independence of matrix elements. 
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In the classical case, the unperturbed energy Eq is not 
constant along a classical chaotic trajectory of the full 
Hamiltonian with a given total energy H — E. Instead, 
it sweeps a range of values, or "energy shell" , Ai?o = Al^, 
and is distributed inside this shell according to a measure 
We{Eq). The form of We{Eq) depends on the form of 
the perturbation V] we will call this measure "ergodic" 
because it is determined by the ergodic (microcanonical) 
measure on the given energy surface H = E. The quan- 
tum analog of this measure characterizes the distribution 
of the "ergodic" eigenfunction (EF) in the unperturbed 
basis. 

Conversely, if we keep the unperturbed energy Eq 
fixed, the bundle of trajectories of the total Hamiltonian 
H, which reach the surface Hq = Eq, has a distribution 
in the total energy E which is described by a measure 
WEoiE). In the quantum case, this measure corresponds 
to the energy spectrum of the Green function (GFS) at 
energy Eq, and has received different names, such as 
"strength function", "local spectral density of states", 
"spectral measure" of the unperturbed eigenstate at en- 
ergy £'0. 

An expression for the latter measure has been given by 
Wigner For a typical perturbation, represented by a 
WBRM, the average measure w{E) = {weo{E)) depends 
on the Wigner parameter , 



(4) 



and has the following limiting forms (see also Refs. 



^ ^El - E\ \E\ < E,,, q » 1 



w{E) 



(5) 



Outside the specified energy intervals, both distributions 
have exponentially small tails. Formulae are valid 
provided pv > I, which is the condition for strong cou- 
pling of neighboring unperturbed states by the pertur- 
bation. In the opposite case pv < 1 the effect of the 
perturbation is small, and we have the so-called pertur- 
bative localization. 

In the limit q >> 1 we have the semicircle (SC) law 
and the width of the energy shell AE — 2Esc — ivV2b — 
4i/2gi?f, ^ Eb where Ei, — h/p is the half width (in 
energy) of the band. In the other limit, q « 1, we 
have the Breit - Wigner (BW) distribution, of width 
AE — 2Ebw — ^Ef,, with the main part inside a width 
F = = 2TTqEb <^ E^. In all these expressions E 

is measured with respect to the center of the distribu- 
tion. Since q << 1 requires pv < Vb, in the BW regime 
the perturbation is not strong enough to couple all states 
within one bandwidth. This means that the BW regime 



corresponds in fact to a sort of partial perturbative lo- 
calization. 

The numerical results presented below are contained 
in the EF matrix Cm„, which connects exact eigenfunc- 
tions "tjjrn , obtained by diagonalization of the Hamiltonian 
matrix (|l|), to unperturbed basis states </>„, 



(6) 



In what follows the eigenvalues Em are ordered, so that 
Em ~ m/p. 

From the matrix Cmn we have found both the statisti- 
cal distribution Wm{n) — Cmn of the eigenstates 4'm on 
the unperturbed ones , and the distribution Wn {m) of 
the unperturbed eigenstates on the exact ones; the mean- 
ing of these distributions is similar to that of the classical 
W and w discussed above. We have then analyzed both 
distributions, and have compared their structures to each 
other and to the SC distribution, paying special atten- 
tion to localization. By localization we shall here mean a 
situation, in which cigenfunctions are localized on a scale 
which is significantly smaller than the maximum one con- 
sistent with energy conservation. Indeed, the size of the 
region which is populated by an eigenfunction (termed lo- 
calization length in the following) is bounded from above 
by the ergodic localization length S'^^ = cpAE, which 
measures the maximum number of basis states coupled 
by the perturbation. This length characterizes the full 
width of the energy shell AE. The factor c depends on 
the definition of localization width (see eqn. (0) below). 
In other words, in a conservative quantum system there 
is always localization in energy, due to the existence of a 
finite AE This fact, which is sometimes a source of 
confusion , is just a trivial consequence of energy conser- 
vation. Here we are interested in localization inside the 
shell 1^, which can be caused by quantum effects. In this 
connection, the matrix size N is an irrelevant parameter, 
provided N ^ d^*^^ is large enough to avoid boundary 
effects. The quantum model (|^) is thus defined by the 3 
physical parameters p, f , and b. 

The localization length d„i of a distribution Wm (n) can 
be defined in several ways. We have used the so-called 
inverse participation ratio (see, e.g., Ref. [^): 
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An) 



(7) 



and similarly for Wn{m). The numerical factor 1/3 ac- 
counts for fluctuations in individual distributions. These 
distributions are assumed to be Gaussian and indepen- 
dent II . 

In order to suppress large fluctuations in individual 
distributions of both types Wm (n) and w„ (to) , we have 
taken averages over 300 of them, chosen around the cen- 
ter of the spectrum. Since different distributions cover 
different regions of the n (respectively, to) space, prior to 
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averaging they have to be shifted into a common region. 
This we have done in two different ways, namely, either 
by counting the site label n in Wmin) starting from the 
center of the energy shell i.e., from the reference site m 
(and vice- versa in the case of w„(m)) (circles in Figs. 1,2), 
or from the center nc(m) of Wm(n), defined by 



nc(m) = ^ Wm{n) 



(8) 



The two types of average will be denoted by 
{W{n)),W{n) respectively. In particular, W{n) yields 
the average shape of an eigenstate (full Hue in Figs. 1,2). 

First we shall discuss the distributions Wm{n)- In ||] 
it was shown that the average localization length d = 
((X^n ^mi^))~^) obeys a scahng law of the form 



1 - 



< 1 



where 



A 



(9) 



(10) 



Here a « 1.2; factor c can be directly calculated from the 
limiting expression (||) for w, which gives c « 0.92. 

The empirical relation (||) has been found |^ to hold 
in the whole interval A < 2.5 and was confirmed in the 
present studies up to A w 7 

The parameter A has been shown j|] to play the role of 
an ergodicity parameter because, when it is large, the 
localization length approaches its maximal value d^'^\ 
which means that the eigenfunctions become ergodic, 
i.e., delocalized over the whole energy shell. Notice 
that in the BW region the ergodic localization length 
d^"') = npT = 27r2&g, and A « ab/2qn'^ > 1 H since 
q ^ 1 (and 6 ^ 1 in quasiclassical situations) . Hence, 
localization is only possible in the parameter range in 
which the local density of states follows the SC law. This 
domain is the main object of the present studies. 

In the case A >> 1 (Fig. la) we have found that the 
averaged distributions {W{n)),W{n) are fairly close to 
the SC law: a remarkable result, because that law was 
theoretically predicted for the other distribution, namely, 
for the GFS spectrum {w{m)). We presume that the de- 
viations from the SC law which are observed in the dis- 
tribution W{n) are due to the not very large value of the 
ergodicity parameter (A — 3.7). The numerical values of 
the localization parameter (j9|) are /3 = 0.94 and (3 = 1.08 
for the two types of average, respectively, in a reasonable 
agreement with the average fid = 0.97 computed from (|^) 
for A = 3.6. For finite q the average distributions of both 
types are bordered by two symmetric steep tails, which 
apparently fall down even faster than the simple expo- 
nential. The structure of these tails will be discussed in 
detail elsewhere. 



The structure of EFs is completely different in the case 
A << 1. (Fig. lb). Whereas individual eigenstates ex- 
hibit large fluctuations, the main part of the average dis- 
tribution ( with respect to the center ric) W{n) shows a 
clear evidence for exponential localization, with localiza- 
tion length in agreement with the empirical formula (|^). 
The width of the main part is small {(3 — 0.24), which 
is again close to average Pd ~ 0.21 for A — 0.24. We 
have found that the main part of the distribution can be 
represented reasonably well by a simple expression: 



W{n) 



2/ Til 



cosh {2n/l) 



(11) 



where the parameter / is related to the localization length 
by / ~ An^'^d. If, instead of averaging the EFs with re- 
spect to their centers, we average them with respect to 
the center of the energy shell, a nice SC (with some tails) 
reappears (Fig. lb, (3 — 0.99) in spite of localization . 
This shows that, in the average, the EFs homogeneously 
fill up the whole energy shell; in other words, their centers 
are randomly scattered whithin the shell (see also fig. 3). 
The latter type of averaging provides a new method for 
calculating ergodic d^*^', and hence the important local- 
ization parameters Pd and A (^). 

Now we turn to the analysis, in the case A <C 1, of 
the other type of distribution: the GFS, or local spec- 
tral density of states w„(m), which is obtained from the 
columns of the matrix Cmn ■ The structure of this distri- 
bution is quite different from that of EFs (represented by 
matrix rows) . Averaging with respect to their centers or 
with respect to the shell center now yields similar results, 
which well fit the SC distribution in both cases (Fig. 2: 
(3 = 0.97 and 0.99, respectively, cf. Fig.lb with (3 = 0.24 
and 0.99). So, GFS look extended, yet they are local- 
ized! This is clear from the average of the corresponding 
individual /3-values: < (3 > — 0.20. The explanation of 
this apparent paradox is that, though each GFS is ex- 
tended over the shell, it is sparse, that is, it contains 
many 'holes'. 

The difference in the structure between EFs and GFS is 
clear from Fig. 3, where solid vertical bars show the main 
parts of EFs. GFS are represented by horizontal dashed 
lines whose sparsity immediately follows from scattered 
localized EFs. Our physical interpretation of the above 
described structure is the following. Spectral sparsity de- 
creases the level density of the operative EFs (that is, the 
ones which are actually excited in a given initial state). 
This is the essential mechanism of quantum localization, 
via decreasing the relaxation time scale Yet, the 

initial diffusion and relaxation are still classical, similar 
to the ergodic case, which requires extended GFS. On 
the other hand, EFs are directly related to the steady- 
state distribution, both being solid because of the homo- 
geneous diffusion during the statistical relaxation. 

In conclusion, we have analyzed the structure of the 
GFS (local spectral density of states) and of the eigen- 



3 



functions for a class of Random Matrices which comes 
much closer to the structure of the real Hamiltonian ma- 
trix of a conservative system than the conventional full 
Random Matrices. Wc have provided numerical evidence 
for the existence of both a delocalized regime, in which 
eigenfunctions have maximal size, with an average shape 
close to the semicircle law, and of a localized regime, in 
which the size of EFs is much smaller than the semicircle 
width. More precisely, quantum localization introduces a 
symmetry breaking, in the sense that the eigenfunctions 
are solid, narrow and randomly scattered inside the en- 
ergy shell, while the GFS remain extended over the whole 
shell but become sparse. In classical language, the lat- 
ter situation means that, although classical trajectories 
are ergodically distributed over the whole energy surface, 
the quantum eigenfunctions cover but a small region of 
the latter. Thus our results indicate that quantum lo- 
calization is a more general phenomenon than commonly 
believed, and suggest similar investigations for realistic 
Hamiltonian, conservative, classically chaotic systems. 
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FIG. 1. Structure of ergodic (la) and localized (lb) eigen- 
functions. Each figure corresponds to a single matrix with 
parameters N = 2560, = 0.1, b = 16,p = 40,g = {pvf/b = 1 
(a), and iV = 2400,u = 0.1,& = 10, p = 300, g = 90 (b). 
The fat full line is the semicircle law (5). Solid lines were ob- 
tained by averaging 300 eigenfunctions with respect to their 
centers; circles, by averaging the same eigenfunctions with 
respect to the centers of their energy shells. In the ergodic 
case (a) A = 3.7, all distributions are close to one another 
apart from fluctuations. In the localized case (b) A = 0.24, 
the average with respect to centers nc{m) of the distributions 
Wm(n) shows a clear localization with {3 = 0.24, while the 
other average remains close to semicircle, with /3 = 0.99. 

FIG. 2. Structure of the GFS (local DOS) for a single ma- 
trix, with the same parameters as for fig. lb. The same aver- 
ages as in Fig.l are shown, and unlike that case they are close 
to each other and to the semicircle law. 

FIG. 3. A comparison of the structure of eigenfunctions 
and of GFS in the localized case of fig. lb. Solid vertical bars 
represent the widths An of individual eigenfunctions over the 
unperturbed basis. Horizontal dotted lines show the size Am 
of the local spectrum for individual basis states. Although 
all basis states have comparable sizes, close to the size of 
the energy shell, they are very sparse {l3 = 0.20), due to the 
fact that EFs are strongly localized, and irregularly scattered 
inside the energy shell. 
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